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Abstract. We report the results of the observation of the nearby satellite galaxy Segue 1 
performed by the MAGIC-I ground-based gamma-ray telescope between November 2008 and 
March 2009 for a total of 43.2 hours. No significant gamma-ray emission was found above the 
background. Differential upper limits on the gamma-ray flux are derived assuming various 
power-law slopes for the possible emission spectrum. Integral upper limits are also calculated 
for several power-law spectra and for different energy thresholds. The values are of the order 
of IQ-^^ ph cm-2 s-i above 100 GeV and 10~i2 

cm s above 200 GeV. Segue 1 is 
currently considered one of the most interesting targets for indirect dark matter searches. In 
these terms, the upper limits have been also interpreted in the context of annihilating dark 
matter particles. For such purpose, we performed a grid scan over a reasonable portion of the 
parameter space for the minimal SuperGravity model and computed the flux upper limit for 
each point separately, taking fully into account the peculiar spectral features of each model. 
We found that in order to match the experimental upper limits with the model predictions, 
a minimum flux boost of 10'^ is required, and that the upper limits are quite dependent on 
the shape of the gamma-ray energy spectrum predicted by each specific model. Finally we 
compared the upper limits with the predictions of some dark matter models able to explain 
the PAMELA rise in the positron ratio, finding that Segue 1 data are in tension with the 
dark matter explanation of the PAMELA spectrum in the case of a dark matter candidate 
annihilating into r^r~ . A complete exclusion however is not possible due to the uncertainties 
in the Segue 1 astrophysical factor. 

Keywords: MAGIC, Segue 1, Dwarf Spheroidal Galaxies, Dark Matter 
ArXiv ePrint: 1103.0477 

1 Introduction 

A major open question for modern physics is the nature of dark matter (DM): strong ex- 
perimental evidences suggest the presence of this elusive component in the energy budget 
of the Universe [see, e.g., 1, 2], without, however, being able to provide conclusive results 
about its nature. One of the most popular scenarios is that of weakly interacting massive 
particles (WIMPs), that includes a large class of non-baryonic DM candidates with a mass 
typically between few tens of GeV and few TeV and an annihilation cross section of the order 
of weak interactions [3, 4]. Natural WIMP candidates are found, e.g., in SUperSYmmetric 
(SUSY) extensions of the Standard Model (SM) [5, 6], and in theories with Universal Extra- 
Dimensions [7-9]. Recently, a number of experimental results have appeared, which may be 
interpreted in terms of DM [10]: in particular, the measurement of the positron fraction in 
cosmic-rays in the 10 — 100 GeV range by the PAMELA satellite [11], and the energy spectra 
of electrons and positrons above 300 GeV measured by Fermi/LAT, H.E.S.S., and ATIC 
[12-14]. Lately, also the search for DM with the use of underground recoil experiments has 
produced interesting results [15-18]. In this paper we will discuss the possibility of detecting 
signatures of DM annihilation in the very high energy band (i.e. above 100 GeV) of the elec- 
tromagnetic spectrum using as particularly convenient target the Segue 1 dwarf Spheroidal 
galaxy (dSph). 
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Assuming the ACDM cosmological model, DM structures form by hierachical collapse 
of small overdensities and are supposed to extend in mass down to the scale of the Earth or 
even below [19-21]. DM structures may also host smaller satellite structures and it has been 
proposed that the dSphs in the Milky Way (MW) may have formed within some of these 
subhalos hosted in the larger MW DM halo [22-25]. So far, around two dozen dSphs have 
been identified in the MW. They represent excellent targets for indirect DM searches due 
to their very large mass-to-light ratios, low baryonic content (which disfavours gamma-ray 
emission from conventional astrophysical sources) and their location often at high galactic 
latitudes, where contamination by Galactic background is subdominant [26-29]. 

To this class belongs Segue 1, a satellite galaxy discovered by Belokurov et al. [30] 
in the Sloan Extension for Galactic Understanding and Exploration [SEGUE, 31]. It is 
located at a distance of 23 it 2 kpc from the Sun (28 kpc from the Galactic Center) at (RA, 
Dec) = (10. 12'^, 16.08°). Despite the fact that its nature has been debated after discovery 
[30, 32-34], it has now been interpreted more clearly as an ultrafaint MW satellite galaxy in 
Ref. [35] through the identification of several new member stars (66 instead of the previous 
24). This interpretation has also been confirmed in Ref. [36] using a Bayesian analysis of star 
membership probabilities. Segue 1 was highlighted as a good target for indirect DM detection 
in Refs. [37, 38]. We will show later on in section 4 that, following our calculations. Segue 1 
results to be the most DM dominated dwarf in our galaxy, known so far. The Fermi/LAT 
[39] has observed Segue 1 in its survey observation mode: results from the first 3 months of 
data have been presented in Refs. [40], while results from the first 9 months of data have been 
presented in Refs. [41, 42]. No gamma-ray signal was detected. Implications for DM models 
have been considered in Refs. [43, 44]. Segue 1 was also observed by the X-ray telescope 
onboard the SWIFT satellite. Data from this short observation (~ 5 ks) have been analyzed 
[45] where no detection is reported either. 

We present here the results of the observation of the Segue 1 satellite performed by the 
Major Atmospheric Gamma-ray Imaging Cherenkov (MAGIC) telescope. The observation 
took place between November 2008 and March 2009, for a total of 43.2 hours (29.4 hours in 
good observational conditions after quality cuts). No significant detection above the back- 
ground was found for energies above 100 GeV. In this paper we describe the data analysis 
and determine upper limits (ULs) for the gamma-ray emission. We extensively study the 
dependence of the ULs on a set of parameters like the expected energy spectrum of the 
source or on experimental features like the energy threshold and the angular resolution. The 
MAGIC-I observation is then used to put constraints on some models of DM. We focus on 
the case of annihilating SUSY DM candidates and compute a grid scan over a reasonable 
range of the SUSY parameters. Point per point we use the predicted energy spectrum to 
compute specific ULs, showing that the dependence on the energy spectrum is even more 
prominent than in the case of power-law spectra. ULs are also given for DM models that fit 
the PAMELA and Fermi/LAT cosmic-ray spectra. 

The paper is organized as follows: in section 2 we describe the observation of Segue 1 
with MAGIC-I and discuss the analysis. In section 3 the calculation of the ULs on the 
differential and integral flux is outlined and results are reported in the case of power-law 
spectra. We focus on a possible DM interpretation in section 4. Conclusions are drawn in 
section 5 together with a brief summary. 
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2 Observation of Segue 1 with the MAGIC-I telescope 



The MAGIC experiment for ground-based gamma-ray astronomy consists of a system of two 
telescopes operating in stereoscopic mode since fall 2009 at the Canary Island of La Palma 
(28.8° N, 17.8° W, 2200 m a.s.l.). However, the observation of Segue 1 was carried out when 
only the first telescope was operating. The data were taken during dark nights between 
November 2008 and March 2009, for a total exposure time of 43.2 hours. The source was 
surveyed at zenith angles between 12.7° and 33.9°, which guarantees a low energy threshold, 
in false-source tracking (wobble) mode [46], in which the pointing direction alternates every 
20 minutes between two positions, offset by ±0.4° in RA from the source. 

The data analysis was performed using the standard MAGIC-I analysis and reconstruc- 
tion software. After the calibration and image cleaning [47, 48], the hadronic background 
rejection is achieved through a multivariate method called Random Forest (RF) [48, 49]. 
The algorithm uses the image Hillas parameters [50] and timing variables to compute a 
gamma-ray/hadron discriminator called hadronness by comparison with Monte Carlo (MC) 
gamma-ray simulations. Hadronness ranges from (for showers confidently identified as ini- 
tiated by gamma-rays) to 1 (for those clearly showing the features of a hadronic cosmic-ray 
initiated shower). The RF method is also used to estimate the energy of the showers. The 
detector MC simulation was tuned to fit the actual telescope performance at the time of the 
observations. After standard quality cuts in trigger rate, atmospheric condition, and a check 
of the distribution of basic image parameters, around 30% of the Segue 1 data were rejected, 
resulting in 29.4 hours of data. The same quality cuts were applied to a sample of Crab 
Nebula data, taken in similar experimental conditions, to cross-check the RF routines and 
the analysis cuts. 

With respect to the standard MAGIC-I analysis, an additional cut was introduced to ac- 
count for the presence of the star r/— Leonis, with apparent magnitude V = 3.5, B — V = 0.02 
and (RA, Dec) = (10. 12^^, 16.76°), located at 0.68° from the position of Segue 1. The light of 
the star entered the trigger region of the MAGIC-I camera and, as a consequence, produced 
local inhomogeneities, which are not treated automatically by the standard analysis at ener- 
gies below 200 GeV. Therefore, all events below this energy and spatially coincident with an 
optimized region around the position of the star in the camera were discarded. With such 
additional star-cuts, around 20% of the events were rejected, but the local inhomogeneities 
caused by the star in the trigger area were efficiently washed out. To check the effect of the 
cuts in the effective area and eventually in the sensitivity, the same star-cuts were applied to 
the Crab Nebula data, and to the MC gamma-ray events for consistency. The final energy 
threshold was at 100 GeV and the star-cuts had the main effect of slightly reducing the 
effective area at lower energies. 

The number of candidates for gamma-ray events from the direction of the source was 
estimated using the distribution of \a\ angles^, which is related to the orientation of the 
showers. The best set of cuts for the hadronness and |a| parameters were optimized on 
the Crab Nebula data sample, separately in logarithmic energy bins (dynamic cuts), and 
under the assumption of a point-like emission. In section 4 we will further discuss this 
issue. In figure 1, the |a|-plot above 100 GeV is shown. For this plot, the number of excess 

^The |a| parameter is defined as tlie angle between the major axis of an image and tfie direction from tlie 
image center of gravity to a reference point in the field-of-view, which is the nominal source position for the 
ON data, and a control background point for the OFF data [50] 
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Figure 1. |q;|— plot from 29.4 hours of observation of Segue 1 above 100 GeV with the MAGIC-I 
telescope. The vertical dashed line represents the fiducial region for the static cut of |a| < 14°, where 
the signal is expected. Red points represent the signal (ON sample), black points the background 
(OFF sample), and green points their difference. 



events was computed in a fixed fiducial signal region (static cut) with \a\ < 14°^ and is 
Nexci> 100 GeV)= —279 ± 329, corresponding to a significance of —0.85 a, computed using 
eq. (17) of Li & Ma [51]. 

In figure 2, we show the significance map for the sky region around Segue 1. For this 
plot the source independent DISP method was used [52]. The energy threshold in this case 
was raised to 200 GeV to cope with the presence of the star in the trigger region. The 
significance distribution is consistent with background fluctuations. 



3 Upper Limits for Power-Law Spectra 

In the previous section we presented the Segue 1 data collected by the MAGIC-I telescope, 
showing that, above 100 GeV, results are consistent with no signal over the background. We 
derive now differential and integral flux ULs for the gamma-ray emission from the source. In 
the computation of the ULs, we will distinguish between "true" energy E (defined for the 
MC gamma-ray simulated events) and "reconstructed" energy Ej-ec, calculated by means of 
the RF method (see section 2) for the real as well as for the MC gamma-ray data. 

^ The static |a|-cut at 14° is used for display purpose. In the computation of the flux upper limits, the 
|a|-cuts are instead optimized more suitably in a dynamic way in each energy bin. In the first two bins, were 
most of the events are concentrated, the optimized |Q!|-cut is 14°. 
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Figure 2. Significance map for events above 200 GeV in the Segue 1 sky region. Tlie black cross 
marks the position of the central core of Segue 1 and the telescope PSF is also shown. The significance 
distribution is consistent with background fluctuations. 



3.1 Differential Upper Limits for power-law spectra 

In order to calculate the differential flux UL in an energy bin AErec one can estimate the 
UL in the number of excess events N^J^ {AErec) following the so-called "Rolke method" [53, 
assuming 30% systematic uncertainty in the signal efficiency and 95% confidence level] from 
the |a|— plot (obtained after the usual analysis cuts and only for reconstructed energies within 
the bin AErec)- The average effective area Aeff{E; AErec) is calculated from the effective 
area for gamma-rays of true energy E, after all analysis cuts, including Erec G AErec, and 
convolved with the energy spectrum. The latter, in this case, is assumed to be a power-law 
of the form S{E) oc (E/E^)^, T being the spectral index and E^ the pivot energy for the 
particular energy bin, defined as: 



J^ESjE) Aeff{E-AErec)dE 

S{E) Aeff{E-AErec)dE ■ 
Finally, the differential flux UL reads as follows: 



d^ 



UL 



1 



(3.1) 



(3.2) 



dE ^ teff Aeff {E; AErec) S{E)dE 

and is computed as the UL in the number of excess events divided by the effective observation 
time and the average effective area weighted with the expected gamma-ray spectrum. Thus, 
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Figure 3. Differential flux upper limits from Segue 1 as in table 1. As reference, the Crab Nebula 
differential flux (red dashed line) [56] and its 10% (blue dashed line) and 1% (yellow dashed line) 
fractions are also drawn. 



results depend on the assumed source spectral shape (although this dependence gets less 
significant as the binning goes finer). The differential UL is measured in ph cm~^ TeV~^. 

Table 1 summarizes the ULs in four reconstructed-energy logarithmic bins between 
100 GeV and 10 TeV using eq. (3.2) and assuming different power-law spectra, with spectral 
indexes T = —1.0, —1.5, —1.8, —2.0, —2.2, —2.4 respectively, as done in Ref. [54]. Moreover, 
the case of F = —1.5 is considered, as a reference for hadronic annihilation DM models [55]. 
The results are also shown in figure 3. 



Table 1. Differential Segue 1 flux upper limits for several power-law gamma-ray spectra in four 
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3.2 Integral Upper Limits for power-law spectra 

Extending the computation of the ULs to energies above a given threshold Eq, we end up 
with the foUowing expression: 

^ " teff J,°^A,ff{E;Erec>Eo)SiE)dE 

where S{E) is again the assumed energy spectrum of the source, Nj/^^{Erec > Eq) is derived, 
as mentioned above, from the corresponding |a|— plot for Erec > Eq, and the effective area 
is reweighted for the assumed source spectrum. 

Table 2 and figure 4 present the integral ULs achieved by the MAGIC-I observation of 
Segue 1 for different energy thresholds and different power-law spectra with spectral index 
r = —1.0, —1.5, —1.8, —2.0, —2.2, —2.4 respectively, as in the previous section. The results 
are comparable with the results of other lACTs on observations of dSphs [57-62]. 



Table 2. Integral Segue 1 flux upper limits for several power-law spectra and different energy thresh- 
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In the case of integral ULs, we see that the dependence on the different power-law 
gamma-ray spectra is slightly stronger than in the differential ULs case, and the UL for the 
r = —1.0 spectrum is at most a factor of two stronger than for the F = —2.4 spectrum for 
low energy thresholds. 

In addition, we clearly see the effect of the statistical fluctuations (characterized by 
the significance of detection (JLi,Ma) on the integral ULs value. The significance of detection 
depends on analysis cuts and the level of residual background: this has quite a strong impact 
on the UL of the number of events N]JJ^, calculated with the Rolke method, and eventually 
on the flux UL (see tables 1 and 2). Quantitatively, a smaller significance corresponds to a 
lower (more stringent) UL: e.g., going from o"Li,Ma = —1 to 1, the UL gets worse by a factor of 
3. This is an intrinsic feature of the stastical methods exploited in the analysis and it should 
be taken into account when comparing ULs from different analyses. 

For this reason, in order to estimate the effect of the spectral slope and of the energy 
threshold on the value of the integral UL, without being biased by the fluctuation due to 
different values of the signiflcance, we computed again the ULs assuming a null value for 
cTLi^Ma^ in particular dashed lines in figure 4 are obtained with a significance equal to zero 
(i.e. with number of ON events equal to the number of OFF events in the signal region of 
the |a|-plot) for different values of F and £"0. 
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Figure 4. Integral flux upper limits from Segue 1. The arrows indicate the integral flux upper limits 
as in table 2 for different power-law speetra and energy thresholds. On the contrary, the dashed lines 
indicate the corresponding integral upper limits if zero significance cTLi^Ma is assumed. 



4 Constraints on Dark Matter Models 

In this section Segue 1 is treated as target for indirect DM searches. Assuming a particular 
form for Segue 1 DM halo, we translate the ULs calculation described before into constraints 
on the DM annihilation rate. 

The gamma-ray flux due to DM annihilations depends on i) the intrinsic DM density 
distribution in the source, ii) the particle physics characteristics of the DM candidate and 
iii) the telescope energy resolution e, the field of view AQ within which the putative signal 
is integrated, and the energy threshold Eq. It is usually factorized in two terms: 

^> Eo,An) = '^^P{> Eo) J{An), (4.1) 

where ^ is the so-called particle physics factor and reads as follows: 

where {cr^nnv) is the velocity averaged annihilation cross section, rriy. is the DM particle mass, 
and "^^=1 dNl^ / dE = dN^/dE is the sum over all the n possible annihilation channels 
producing photons is the particular branching ratio for channel i). 
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The term J(Ai}) is the astrophysical factor and is given by the hne-of-sight integral 
over the DM density squared within a sohd angle Afi: 

J{An)= [ [ p'^{r{s,n))dsdn. (4.3) 
J An J los 

Motivated by results from cosmological simulations, to describe the Segue 1 DM density 
distribution we used the Einasto profile [63]: 

mN(r)=p.e-2-N-)^'"-^], (4.4) 

which produced a good fit to the subhalos simulated by the most recent A^-body simulations 
[21]. Eq. (4.4) is defined by three parameters: the scale density ps, the scale radius and 
the index n, which typically ranges from —3 to 7. For the computation of the astrophysical 
factor, we have used ps = 1-1 x 10*^ MQkpc"'^, = 0.15 kpc and n = 3.3, chosen among the 
central values after the marginalization of the likelihood used in Ref . [44] . The astrophysical 
factor is compatible with the recent estimation of Ref. [64]. The astrophysical uncertainty 
on J(Ar2) has been estimated to be slightly larger than one order of magnitude at 2(T level 
[44]. 

In figure 5 we show how the gamma-ray fiux increases when integrating within larger 
regions around the center of the source. We see that the source is slightly extended compared 
to a representative MAGIC-I angular resolution of 0.1°. On the other hand, the contamina- 
tion of the signal in the background region, calculated at a distance of 0.8° from the "signal" 
region and with the same extension, is well below 1%. The energy dependent \a\ cuts applied 
in the analysis determine the angle of integration which fix the fraction of the total astrophys- 
ical factor that must be considered for the ULs computation. In order to get this angle we 
have produced a toy MC simulation for extended sources starting from point-like MC events 
and spreading the original gamma-ray arrival directions according to the source DM density 
distribution of figure 5. Since the source is only slightly extended and the assumed luminosity 
is very peaked at the center, this method mimics reasonably well a dedicated analysis for an 
extended source. We have found that the applied \a\ cuts correspond to an angular integra- 
tion of 0.14° from the Segue 1 center (a solid angle of 10~^ sr). From figure 5, one can see 
that such angle of integration encloses 64% of the total astrophysical factor. Therefore, with 
the values mentioned above for scale radius and density which leads to a total astrophysical 
factor of J(AJ7) = 1.78 x 10^^ GeV^ cm~^ sr, the effective astrophysical factor within the 
considered analysis cuts is J(Af]) = 1.14 x 10^^ GeV^ cm~^ sr. This is, to our knowledge, 
the largest astrophysical factor among the known dwarf galaxies. We compared in fact our 
results with the astrophysical factors of table 4 of Ref. [54], taking care of integrating over 
the same region. For a solid angle of 2.4 x 10~^ sr, we get 95% of the total astrophysical 
factor of Segue 1, i.e. J{AQ) = 1.69 x 10^^ GeV^ cm~^ sr, while the largest value in Ref. [54] 
is 1.2 X 10^^ GeV^ cm~^ sr in the case of Draco. 

Focusing now on the particle physics factor, we restrict ourselves for now to the case of 
a SUSY model in which the presence of a discrete symmetry (i?— parity) guarantees that the 
Lightest Super symmetric Particle (LSP) is stable over cosmological timescales and, therefore, 
a good WIMP candidate. We will consider a 5-dimensional subspace of the Minimal Super- 
symmetric Standard Model (MSSM) called mSUGRA [65, 66], which is defined by universal 
masses for the gauginos (?7ii/2)) scalars (mg), and trilinear couplings (^o)) as well as by the 
ratio of the vacuum expectation values of the two Higgs fields (tan /3) and the sign of the 
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Figure 5. Fraction for squared DM density around Segue 1 integrated within an angular opening 
with respect to the total squared DM density, as a function of 'I' itself. The dashed line indicates 
the angular opening corresponding to the convolution of the energy dependent |a|— cuts used in the 
current analysis, which were optimized for point-like sources. An Einasto profile, as defined in the 
text, is assumed to describe the DM density distribution. 



Higgsino mass term sign(/x). In the majority of mSUGRA models, the LSP is the lightest 
neutralino Xi ^ linear combination of the super-partners of the gauge bosons and neutral 
Higgs bosons. 

In order to study the phenomenology of mSUGRA we performed a grid scan over the 
parameter space. We spanned (with linear steps) the regions indicated in table 3, chosen 
to provide neutralino masses within the range of detection for MAGIC-I. For each direction 
we considered 40 steps and in each bin along mo and mi/2 we randomly selected a point 
within that bin. The scan has been done once with a positive sign(^) and then again with a 
negative sign(;u), for a total of 5 x 10^ points. For each mSUGRA model we used DarkSUSY 
5.0.4 [67] (which includes the virtual internal Bremmstrahlung effect [68]) i) to test if the 
model is physical, ii) to check if it passes the Standard Model (SM) experimental constraints 
implemented in the code (e.g. LEP bounds on Higgs mass nih > 114 GeV, on chargino mass 

> 103.5 GeV and constraints from b sj), and in) to compute the relic density [with 
Isasugra 7.78 69]. 

All the models in the scan that correspond to a neutralino with a relic density compatible 
with the value derived by WMAP data within three times its experimental error cJwmap [1] 
are plotted as crosses in figure 6. The crosses approximatively cover the mass range between 
100 GeV and 1 TeV, due to the range that we considered for tjiq and mi/2 in the scan. On the 



- 10 - 



Table 3. Parameter space for the scan over mSUGRA. From Komatsu et al. [1] Q'^^^^h'^ = 0.1123 
and (TwMAP = 0.0035. 



Parameter 


Range 


Steps 


mo 


OU , OUUU KjQ V 


/in 


mi/2 


; 5000 GeV 


40 


tan/3 


2; 62 


40 




-7000 ; 7000 


40 


sign(^) 


+;- 


2 


Total numb 


er of models scanned 


5.12 X 10*^ 


<^ 


c passing SM bounds 


2.42 X 10*^ 


Sz with r^DM^^ - 


QWMAPL2 ^ q 


42427 


& with r^DM^^ — 


S'dM I ^ "JOWMAP 


4180 



contrary, the values for the cross section span many orders of magnitude. They are mainly 
concentrated around a value of 10^^^ cm^ s~^, but there are models that, due to particular 
mechanisms, are characterized by lower values for the cross section. One example is the 
prominent "strip" that crosses the plane from 10~^^ cm^ to around 3 x 10~^^ cm'^ 
for neutralinos that co- annihilate with stops and staus, or the "tail" at low masses (around 
50 GeV). We stress here that the density of the points in figures 5 and 6 has no probabilistic 
meaning and that the region of points compatible with WMAP may extend until few times 
10~^^ cm^ s~^ if the range of the scan is extended [e.g. figure 25 of 70]. 

For each DM model in the scan, the integral flux UL <^^^{> Eq) can be computed 
following eq. (3.3), using the Segue 1 data and the specific gamma-ray spectrum of the indi- 
vidual DM model. To ensure a more direct comparison with the particle physics predictions, 
we present our results in terms of ULs on the averaged cross section {(t^^^v)^^ ^ following 
eq. (4.2): 



^jj^ _ 87rm^$^^(> Eq 
Jim '^dE 



In Sec. 3.2 we noticed that the integral ULs may change as a function of the energy 
threshold, with lower (more stringent) ULs if Eq is larger than the experimental one [see also 
72]. While this variation is more predictable in the case of power-laws, the situation may 
be less clear for annihilation spectra that contain features and terminate at the DM mass. 
An example of some gamma-ray spectra used in this analysis for a DM mass of 1 TeV is 
shown in figure 7. Therefore, for each model we computed the UL for different values of 
energy thresholds Eq among those listed in table 2. We start by considering the case of null 
significance fiLi^Ma = 0, not to be biased by fiuctuations in the real (TLi,Ma defined above. For 
each annihilation channel and DM mass, we determine the optimal energy threshold Eopt as 
the one with the most stringent flux UL among the set of Eq considered. 

The effect of the energy threshold optimization can be seen in figure 6. In the bottom 
part of the plot, we show the models of the scan compatible with WMAP bounds as black 
crosses. In addition, two representative subsets are also shown using a different color coding 
according to their main annihilation channel (red points for branching ratio B{bb) > 0.85, 
and blue points for B{t~^t~) > 0.7), which are also representatives of a soft and hard gamma- 
ray spectrum respectively (see figure 7). On the top part of the plot, different exclusion lines, 
obtained from the integral ULs, are shown. The blue dashed line is the exclusion curve 
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Figure 6. Annihilation cross section ULs from Segue 1 MAGIC data considering neutralino annihi- 
lating entirely into bb or into r+r". mSUGRA models with a relic density within Scwmap from the 
WMAP value are plotted (black crosses). Among these, neutralinos annihilating mainly in bb and 
T+r" are indicated with red points and blue points respectively. The solid red line indicates ULs for 
a neutralino annihilating entirely into bb while the blue lines the case of annihilations into r+r^. The 
dashed blue line is calculated for a Eo = 100 GeV and by imposing a null significance. The dotted 
blue line is also calculated from a null significance but above the optimized energy threshold Eopt 
as explained the text. Finally, the thick solid blue line represents the UL calculated with the real 
significance and above Eopt. The Egpt is optimized for all the DM masses. Finally, for annihilations 
into T+r", the blue band covers the 2a uncertainty on Je(A17). 

obtained with the integral ULs for the t'^t~ annihilation channel, for (7Li,Ma = and with 
Eq = 100 GeV. The blue dotted line represents the same curve but now recalculated with 
the optimized energy threshold method described above. One can clearly see that the opti- 
mization works at moderate and high DM masses, because at low DM masses, by increasing 
the energy threshold one also loses too many photons. We recall that these two curves do 
not represent real data, because are calculated with (TLi,Ma = 0. The thick blue solid line is fi- 
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Figure 7. Spectra of different annihilation channels. For the channels 66, W^+VF ,r+r we 
used the fits from Ref. [71] . For (pcj) 4e we used the parameterization in Ref. [37] with = 1 GeV. 



nally the exclusion curve calculated with the real significance o"Li,Ma and the optimized energy 
threshold. One can now see that below roughly 1.2 TeV the curve is more constraining than 
the thin blue solid one because in this regime, the best energy threshold is between 100 and 
200 GeV, where the significance is negative (from table 2). Above 1.2 TeV, the optimized 
energy threshold is 251 GeV corresponding to a positive value of o"Li,Ma) and therefore the 
two mentioned lines cross. Finally, the red solid line indicates the UL (with optimized energy 
threshold) in the case of annihilations only into bb. In this case, due to the soft spectrum of 
this channel, the optimized energy threshold is almost all the time 100 GeV, apart from the 
very last points. This is reasonable, because soft spectra do not gain more by selecting higher 
energy photons. In the same figure, it is also possible to see how the ULs depend on the 
shape of the energy spectrum: differences can be larger than one order of magnitude [see also 
54] and, as expected, hard spectra are more constraining. Finally, in the case of annihilations 
into T~^T~ , we also indicate (by means of the blue band) the uncertainty in the astrophysi- 
cal factor at the 2a level, showing how strongly it affects our ULs and prospects for detection. 

In figure 8 we plot again the (cJann'w) predictions (same as in figure 6) for the points 
in the scan (full circles), together with the corresponding ULs calculated above Eopt with 
the method described above (full squares). In addition, all the models in the scan that 
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Figure 8. Annihilation cross section ULs from Segue 1 MAGIC data computed for individual points 
in the scan. Grey crosses indicate the annihilation cross section value for those points in the scan (see 
table 3) that pass the SM constraints and with a relic density lower than WMAP bound. The full 
circles only consider models within Sctwmap from WMAP. For each of these full circles the UL on the 
cross section can be computed from the Segue 1 data (after energy threshold optimization) and it is 
indicated here by a square. Circles and squares are color coded in terms of the enhancement factor 
(see eq. (4.6)). 



correspond to a neutralino with a relic density lower than the value derived by WMAP data 
plus three times its experimental error cJwmap [1] are also plotted as grey crosses. In this 
case, since the standard freeze-out predicts a too small relic density, we would require a 
non-thermal production mechanism for the DM to recover the correct relic density. We are 
not considering the case of multi-component DM: if the neutralino is responsible for only a 
fraction of the total DM density, it would be reasonable to rescale the value of the Segue 1 
astrophysical factor (and consequenty the UL) by the same fraction. It can be seen that, 
for large neutralino masses, ULs concentrate around 10~^^ cm^ s~^, while at lower masses, 
the distribution is wider, because models with a low gamma-ray flux above the threshold are 
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able to produce only loose ULs (as in the case of squares at low masses with UL of the order 
of 10"^'^ — 10~^^ cm^ We stress that each point in the scan should be compared to its 

own UL, and hence the apparent overlap in figure 8 between the ULs (plotted as squares) 
and some models with a relic density below the WMAP value (grey crosses) does not imply 
that any of these models is excluded. In order to avoid this possible misunderstanding we 
decide to compute point per point what we call enhancement factors (ENFs) defined as the 
ratio between the UL on the averaged cross section and the value predicted by mSUGRA: 

ENF = {(Tannvf^/iaannV). (4.6) 

This indicates how much the cross section of the particular model should be increased in 
order to make it detectable. In these terms, one can easily understand which points can be 
excluded by MAGIC-I data on Segue 1 since they would be associated to an ENF smaller 
than one. In figure 8 the color coding is chosen in terms of the ENFs: yellow for points with 
an ENF smaller than 10^, orange for models with 10^ <ENF< 10^ red if 10^ <ENF< 10^ 
and brown if ENF> 10^. 

In figure 9 the ENFs are also plotted as a function of the mass for models compatible 
with WMAP value of the relic density (red crosses) and below (black crosses). The panel in 
the upper right of the figure indicates the ENF distribution for the two sets of models. For 
those compatible with WMAP (red crosses), the lowest ENF is of the order of 10"^. Figure 8 
tells us that the models charactized by the lowest ENF are those with the largest annihilation 
cross section and a neutralino mass above 200 GeV. On the contrary, the majority of the 
points in the scan have an ENF > lO'^. Moreover, it can also be seen that the distribution 
of ENFs is quite wide. As commented before, this large spread is due to the very high (less 
constraining) ULs relative to models with small neutralino masses, as already pointed out 
in Refs. [58, 59], and in general to models with a low gamma-ray flux above the energy 
threshold. In the case of models with a relic density below the WMAP value (black crosses), 
the situation is slightly better: their intrinsic higher cross sections make them closer to their 
ULs, the ENF peaks at values somewhat lower than those for points compatible with WMAP 
and the distribution extends to lower values. We have around 160 points with ENF< 20, 
40 with ENF< 10 and 5 with ENF< 5 even if the exact numbers depend on the number of 
points scanned. 

For the sake of completeness, we point out that the contribution of monochromatic lines 
to the annihilation spectrum is neglected in this study, due to the fact that their contribution 
is subdominant with respect to the continuum emission [73]. 

4.1 Impact on PAMELA preferred region 

Despite the fact that several astrophysical explanations have been proposed to describe the 
rise in the energy spectrum of the positron fraction /{e~^ + e~) measured by PAMELA 
[11, 74-76], an interpretation in terms of DM is still possible [see, e.g., 77-84]. Typically, 
the DM particle has to be heavy and annihilate mainly to leptons. One possibility that has 
been largely studied is that the annihilation to leptons occurs through the production of an 
intermediate state (p, mediator of a new, long range, attractive force [78]. However, a very 
large annihilation cross section is required, about a factor 100 — 1000 larger than the canonical 
value derived for thermal production of {uannv) ~ 10~^^ cm^ s~^. In the following, we test 
our ULs against some of the models proposed in the literature that fit the PAMELA data. 
The regions in the {m-^, {(^annv)) plane that provide a good fit to the PAMELA data are 
shown in figure 10 for a DM candidate annihilating into ;U^/.f~, t~^t~ and for the case of the 
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Figure 9. Enhancement factors as a function of the DM mass for models in the scan providing a 
rcUc density compatible with WMAP value (red crosses) or below (black crosses). The panel in the 
upper right part indicates the distribution of the ENFs for these two sets of models with the same 
color coding. 

intermediate state (j) decaying to e'^e", with = 1 GeV. These regions have been adapted 
from Ref. [83] after rescaling from a local DM density of 0.3 GeV/cm'^ to 0.43 GeV/cm'^ [85]. 

Using again the specific DM annihilation spectra, we plot in figure 10 the ULs obtained 
from the Segue 1 data. We can see that, in this case, the ENFs needed to meet the PAMELA- 
favoured region are much smaller than for mSUGRA, and in the case of annihilation into 
T~^T~ our ULs are probing the relevant regions. However, we recall that the uncertainty 
in the astrophysical factor (see figure 6) is large: if future data on Segue 1 point to an 
astrophysical factor close to the upper end of the currently allowed range, our observations 
of Segue 1 might be able to confirm the exclusion of the PAMELA region for DM particles 
annihilating in r+r", at least for massive DM candidates. Note that for the annihilation into 
— 7> 4e, an additional Sommerfeld enhancement due to the lower DM velocity dispersion 
in Segue 1 compared to the DM velocity dispersion in the MW halo may effectively shift up 
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Figure 10. Exclusion lines for a neutralino DM annihilating exclusively into fJ,'^fJ,~ (green lines) 
or T+T" (blue line) and for a DM candidate interacting with a light intermediate state </> decaying 
into a pair of electrons (pink line). The same annihilation channels (with the same color coding) are 
considered to draw the regions in the plane that provide a good fit to the PAMELA measurement of 
the energy spectrum of the positron fraction. The regions are taken from Ref. [44] , which are adapted 
from Rcf. [83]. We used an astrophysical factor of J(Afi) = 1.14 x 10^^ GcV^ cm~^ sr. 



the PAMELA preferred region by an order of magnitude, at least for low enough <j) masses 
(see also Ref. [44]). 

5 Summary and discussion 

In the present paper we presented the observation of the ultra-faint satellite galaxy Segue 1 
performed by the MAGIC-I telescope (single telescope mode) over 29.4 hours of selected 
data. A description of the data analysis is reported with the result of no detection above the 
background for energies larger than 100 GeV (sections 2). 
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The result is used to compute ULs on the gamma-ray emission from the source, assuming 
different power-law energy spectra. The computation is done first in energy bins (table 1 and 
figure 3) and then for integral ULs above different energy thresholds (table 2 and figure 4). In 
all cases, we averaged the effective area reweighting it with the specific gamma-ray spectrum, 
which allowed us to determine how much the ULs depend on the specific spectrum. We 
also pointed out the fact that one can get more stringent ULs if computed above energies 
larger than the experimental energy threshold, as a result of the interplay between the larger 
sensitivity of the experiment at moderate energies and the assumed spectrum. 

We focused then on indirect detection of DM and produced ULs on the annihilation cross 
section for a large scan of neutralino models within the mSUGRA scenario (figure 8). The 
ULs are derived separately for each point in the scan in order to completely account for the 
dependence on the specific spectral shape. Results indicate that ULs are quite dependent on 
the energy spectrum and a general exclusion plot cannot be drawn to constrain the parameter 
space (figure 8). For this reason, we find it quite useful to provide the results in terms of 
enhancement factors, defined as the intrinsic flux boost needed to meet detection (eq. 4.6). 
Results are shown in figure 9. A mininum boost is found of the order of 10'^ (for models 
compatible with WMAP) while "typical" values are at 10^~^. However, if we loosen the 
constraint and request only that the SUSY models do not overshoot the WMAP value for 
the relic density, then the situation improves since we can have ENFs as low as a few. In 
these terms MAGIC-I data on Segue 1 are not so far from excluding portions of mSUGRA. 

We have also discussed how MAGIC-I data on Segue 1 can be used to test the PAMELA 
results on cosmic-rays. Current results are probing the PAMELA preferred region for the 
case that DM annihilates into r^r". Future improved measurements of the kinematics of 
the stars in Segue 1 may decrease the uncertainty in the line-of-sight integral over the DM 
density squared and rule out this region. 

The robustness of our results depends mainly on the assumptions on the astrophysical 
factor, since an uncertainty of two orders of magnitude (at 2a) remains, as estimated in 
Ref. [44]. Our result can be considered conservative since they do not take into account 
intrinsic contributions to the fiux from the presence of substructures in the "smooth" DM 
halo profile of Segue 1 (which is already "per se" a substructure of the larger MW halo). 
The effects of substructures most likely increase the flux by a factor of a few at the most 
[20, 21, 38, 64]. A second contribution to the flux, unaccounted for in our calculation, comes 
from a particle physics mechanism known as the Sommerfeld effect, which may additionally 
boost up the predictions for the gamma-ray emission from DM annihilation. Its effect on 
the predicted fluxes from halo substructures have been studied in Refs. [44, 86]. The most 
spectacular effects are present for some resonant values of the DM mass where the flux can be 
increased by a factor of ~ 10^. How these models can be constrained by dSphs observation 
has been studied in Refs. [44, 87]. 

Moreover, secondary gamma-ray emission can be expected, that may enhance the pre- 
dicted flux from Segue 1 and make our ULs more constraining. As studied in Ref. [88], the 
main mechanism for secondary emission is Inverse Compton of electrons and positrons pro- 
duced in DM annihilations with CMB and starlight photons. The contribution of the Inverse 
Compton with the CMB has been computed in the case of dSphs in Ref. [54] where it can 
be seen that the UL on the annihilation cross section can improve of one order of magnitude 
for a DM mass around 1 TeV. 

It is beyond the goal of this paper to discuss constraints on other DM scenarios. As 
a brief comment, we simply recall that models with Universal Extra Dimensions predict a 
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stable particle whose annihilation signatures are typically harder than the ones in mSUGRA. 
For these models, one can use the ULs we presented for the t~^t~ channel, which can be taken 
as reference for hard spectra, and therefore draw conclusions similar to ours. Furthermore, 
for the case of decaying DM (that can also fit the PAMELA and Fermi/LAT-HESS results), 
one can rescale our results accounting for a different astrophysical factor (depending on the 
DM density and not on the DM density squared) and constrain the decay timescale instead 
of the annihilation cross section. 

Comparing MAGIC-I results on Segue 1 with those of Fermi/LAT, the different energy 
range covered by the two experiments implies that the latter is more constraining for low 
mass DM candidates, while MAGIC (or Cherenkov telescopes in general) can produce better 
ULs only for DM heavier than few hundreds of GeV. For low mass neutralinos, it is already 
possible for Fermi /L AT to exclude some of the mSUGRA models with a relic density smaller 
than what was measured by WMAP [54]. With 5 years of data, Fermi/LAT may be able to 
probe some of the points of the parameter space considered here. In case of no detection, 
this will have the effect of excluding most of the models with large ENFs, shown in figure 9, 
which normally correspond to low-mass neutralinos. On the other hand, Fermi/LAT will not 
probe the larger DM masses shown in figure 8 and figure 9, where lACTs are more sensitive. 
In this direction, we recall that with the new stereoscopic system, MAGIC-stereo allows for 
improved background rejection specially below 100 GeV (and therefore higher sensitivity), 
improved energy and angular resolutions, and a lower energy threshold [89]. 

In conclusion, while the observations presented here did not result in a detection, and 
the ULs require still quite high (and in some cases unmotivated) fiux enhancement factors 
to actually match the experiment sensitivity, an analysis like the one presented here is able 
to point out details and features that can be important for future deep exposures of this or 
similar objects, using facilities with much improved sensitivity, e.g. the planned Cherenkov 
Telescope Array [90]. 
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